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SUMMARY 


Part  III  of  this  series  of  papers  developed  the  theory  of  high- 
eccentricity  orbits  (e  >0.2)  in  an  atmosphere  having  an  exponential  varia¬ 
tion  of  air  density  with  height,  that  is,  with  the  density  scale  height  H 
taken  as  constant.  Part  IV  derived  the  appropriate  theory  for  low- 
eccentricity  orbits  (e  <0.2)  in  a  more  realistic  atmosphere  where  H 
varies  linearly  with  height  y  (and  (u)  =  dH/dy  <0.2).  The  present  Report 
treats  the  orbits  of  Part  III  when  they  meet  the  air  drag  specified  by  the 
atmospheric  model  of  Part  IV.  Equations  are  derived  showing  how  the  peri¬ 
gee  height  varies  with  eccentricity,  and  the  eccentricity  varies  with  time, 
over  the  major  part  of  the  satellite's  life.  It  is  shown  that  the  theory 


of  Part  III  remains  valid,  to  order 
height  above  perigee.  •' 

j 


,  if  H  is  evaluated  at  a  specific 
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INTRODUCTION 


Part  I  of  this  series  of  papers^  gave  the  theory  for  the  contraction  of  Earth 
satellite  orbits  of  eccentricity  e  <  0.2  under  the  influence  of  air  drag  in  a  spheri¬ 
cally  symmetrical  atmosphere,  over  a  single  revolution  and  also  over  the  complete  life- 

2 

time.  Part  II  did  the  same  for  an  oblate  atmosphere  .  In  Part  III  the  theory  of  Part  I 

3 

was  extended  to  orbits  of  high  eccentricity,  0.2  ^  e  <  1  .  In  Parts  I  to  III  the  air 

density  was  assumed  to  vary  exponentially  with  height,  having  a  constant  density  scale 

4 

height  H  :  in  Part  IV  the  theory  was  developed  ,  for  e  <  0.2  ,  on  the  assumption  that 
instead  of  being  constant,  H  was  a  linear  function  of  height  y  .  Parts  V  and  VI  took 
account  of  the  day-to-night  variation  in  air  density^'^,  again  for  e  <  0.2  ,  and  revert¬ 
ing  to  constant  H  ;  Part  VI  was  devoted  to  the  special  conditions  for  near-circular 

7  8 

orbits.  Subsequently  the  theory  has  been  extended  by  Swinerd  and  Boulton  ’  to  cover  the 
combined  effects  of  oblateness,  day-to-night  variation  and  variation  of  H  with  height, 
for  e  <  0.2  ,  over  a  single  revolution. 

The  present  Report  can  be  looked  on  as  a  marriage  between  Part  III  and  Part  IV,  and 
applies  for  orbits  of  high  eccentricity  (0.2  ^  e  <  1)  in  an  atmosphere  where  H  varies 
linearly  with  height.  The  theory  is  required  for  three  reasons: 

(a)  In  the  real  atmosphere  H  is  not  constant,  but  does  vary  almost  linearly 
with  height;  so  the  constant-H  theory  of  Part  III  may  be  inadequate  in  practice, 
and  the  extent  of  the  inadequacies  is  unknown  until  the  new  theory  has  been 
developed. 

(b)  The  theory  is  needed  in  interpreting  the  measurements  of  geophysical  para¬ 
meters  from  orbital  changes  for  orbits  with  eccentricity  greater  than  0.2:  without 
it,  the  height  at  which  the  measurements  apply  remains  uncertain. 

9 

(c)  In  orbit  determination  programs,  such  as  PROP  ,  a  constant-H  theory  is  used 
in  the  model,  and  the  new  theory  will  allow  the  best  constant  value  of  H  to  be 
selected  for  a  particular  orbit  determination. 

2  ASSUMPTIONS 

The  theory  is  developed  under  the  nine  assumptions  listed  below.  These  assumptions 
are  the  same  as  in  Part  I,  except  for  (c)  and  (h)  and  a  minor  modification  in  (e) . 

(a)  The  atmosphere  is  spherically  symmetrical. 

(b)  The  air  density  p  at  a  given  distance  from  the  Earth's  centre  does  not  vary 
with  time. 

(c)  The  air  density  p  and  the  quantity  -  p/(dp/dr)  ,  which  will  be  called  the 
density  scale  height  and  denoted  by  H  ,  vary  with  distance  r  from  the  Earth's 
centre  in  the  manner  specified  in  section  3. 

(d)  The  resultant  aerodynamic  force  on  the  satellite  acts  in  the  direction 

opposite  to  the  velocity  V  of  the  satellite  relative  to  the  ambient  air,  and  may 
2 

be  taken  as  JpV  SC^  ,  where  is  a  drag  coefficient  based  on  an  effective 

cross-sectional  area  S  .  and  SCL  is  assumed  constant. 


(e)  The  atmosphere  rotates  at  a  constant  angular  velocity  w  not  very  different 
from  that  of  the  Earth. 

(f)  The  Earth's  gravitational  potential  is  taken  as  that  of  a  point  mass  at  its 
centre,  so  that  the  unperturbed  orbit  is  an  ellipse  in  a  plane  passing  through  the 
Earth's  centre. 

(g)  During  one  revolution  the  action  of  air  drag  changes  the  orbit  by  only  a 
small  amount,  whose  square  can  be  neglected. 

(h)  The  orbital  eccentricity  e  lies  within  the  limits  0.2  e  <  1  . 

(i)  Lunisolar  perturbations  are  ignored. 

The  assumptions  have  been  discussed  in  Part  I,  so  little  more  needs  to  be  said 
here,  but  three  points  do  deserve  mention.  First,  high-eccentricity  orbits  are  not  much 
affected  by  atmospheric  oblateness,  so  that  assumption  (a)  is  not  objectionable.  Second, 
the  neglect  of  the  asymmetrical  components  of  the  Earth's  gravitational  field,  particu¬ 
larly  the  third  harmonic,  is  questionable  and  it  is  intended  to  develop  methods  for 
taking  account  of  this  effect  in  a  subsequent  paper.  (For  a  spherically  symmetrical 
atmosphere  the  effects  of  the  much  larger  second  harmonic  are  not  important.)  The  third, 
and  most  important,  of  the  simplifications  is  the  neglect  of  long-term  lunisolar  pertur¬ 
bations,  although  their  effect  on  the  single-revolution  results  is  negligible.  The 
problem  arises  over  months  or  years,  when  lunisolar  perturbations  can  be  very  large  for 
high-eccentricity  orbits.  However,  the  effects  are  then  to  some  extent  self-eliminating, 
because  the  perigee  will  either  rise  to  heights  where  the  drag  is  obviously  negligible, 
or  will  be  forced  down  to  a  height  where  the  satellite  will  decay.  Even  in  the  long 
term,  lunisolar  perturbations  to  the  perigee  height  are  generally  small  for  orbits  of 
eccentricity  less  than  about  0.5,  and  the  theory  is  often  also  still  applicable  at  higher 
eccentricities.  At  eccentricities  higher  than  0.8,  however,  lunisolar  perturbations  are 
generally  dominant  and  we  shall  regard  0.8  as  the  likely  practical  upper  limit  for  e  . 

3  AIR  DENSITY 

The  model  for  the  variation  of  density  p  with  distance  r  from  the  Earth's 
centre  is  as  specified  in  Part  IV.  The  density  scale  height  H  is  defined  by  the 
equation 

"  ■  -jjfe  •  <° 

and  if  H  were  constant  equation  (1)  could  be  integrated  to  give  the  simple  exponential 
variation  of  density  with  height, 

p  =  0p  exp( - ,  (2) 

where  suffix  p  denotes  conditions  at  perigee,  and  in  practice  equation  (2)  is  used  only 
for  r  >  r 

P 


Here  we  assume  instead  that 


p  =  op  1  +  b 


(r '  rp)2J  '-p(-  v*)  ■ 


(3) 


where  H  and  b  are  taken  constant  during  one  revolution.  On  substituting  (3)  into 
P  2 

(1),  and  vrriting  a  =  2bHp  and  ^r  -  r^yH^  =  n  ,  we  find 


H  =  H 


a(r  '  o 

H  +  - ^ 

=  H  +  p(r  -  r  ) 

p  1  -  a(n  -  in  ) 

P  '  P/ 

say, 


(4) 


where 


2bH 


p  1  +  u(n  -  in2) 


(5) 


At  heights  less  than  2.4  Hp  above  perigee,  ip  for  n  <  2.4  ,  the  term  (n  -  $n^)  ^  J  , 
and  if  p  is  small  it  will  usually  be  an  adequate  approximation  to  take  a  =  p  .  (The 
situation  for  n  >  2.4  is  discussed  in  section  5.) 


In  practice  p  is  usually  of  order  0.1:  Fig  1  shows  the  variation  of  H  and  of  p 
given  by10  the  COSPAR  International  Reference  Atmosphere  1372,  for  two  levels  of  solar 
activity.  It  is  seen  that  a  linear  variation  of  H  over  a  height  range  of  50-100  km  is 
a  satisfactory  approximation  and  that  p  does  not  exceed  0.2  unless  the  height  falls 
below  190  km  (low  solar  activity)  o.  230  km  (high  solar  activity).  In  developing  the 
theory  in  this  Report  we  shall  assume  p  <  0.2  . 


It  may  seem  simplest  to  evaluate  p  at  perigee  height.  But  the  subsequent 
analysis  shows  that  the  value  of  H  at  a  height  1.5  Hp  above  perigee  is  of  special 
importance;  so  it  is  recommended  that  the  numerical  value  of  p  at  a  height  about  1.5  H 
above  perigee  should  be  used.  With  this  rule,  Fig  1  shows  that  p  <  0.2  if  the  perigee 
height  exceeds  160  km  (low  solar  activity)  or  190  km  (high  solar  activity). 


P 


The  expression  (3)  for  p  will  be  used  for  integrations  over  one  revolution,  to 
obtain  the  changes  in  the  orbital  elements  over  one  revolution.  These  equations  are  then 
integrated  over  a  long  time-interval,  perhaps  nearly  all  the  satellite's  lifetime. 

During  this  time  there  will  be  a  variation  in  the  perigee  height  and  therefore  in  Hp  . 
For  this  further  integration  we  assume,  in  accordance  with  (4),  that 


H 

P 


’  P° 


+  pfr 
1  P 


(6) 


where  the  second  suffix  0  denotes  initial  values.  If  0.2  e  4  0.8  ,  it  is  found  that 

0  .k  r  „  -  r  0.5  H  .  .  Thus  the  change  in  H  due  to  the  p  term  is  never  more  than 
pQ  P  p0  p 

10%  if  p  <  0.2  ,  and  the  error  of  order  Jp  between  b  and  p  in  equation  (5)  leads 

to  an  error  of  less  than  1%  in  H  .  This  is  acceptable,  because  H  is  rarely  known 

P  P 

with  an  accuracy  better  than  2%.  We  shall  therefore  normally  utilize  equation  (5)  in  the 
2 

form  b  =  p/2H  ,  although  the  fuller  version  is  needed  at  one  point  in  the  analysis. 
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4  BASIC  EQUATIONS 
4. 1  Air  density 

If  a  is  the  semi-major  axis  of  the  orbit,  the  radial  distance  r  of  a  satellite 


may  be  expressed  in  terms  of  the  eccentric  anomaly  E  as 

r  =  a(1  -  e  cos  E)  , 


so  that 


r  -  r 


ae(1  -  cos  E)  , 


(7) 

(8) 


and  (3)  becomes 


p  |l  +  bx2(1  -  cos  E)2 j  exp z(1  -  cos  E)|  , 


where 


x  =  ae 


(9) 

(10) 


and 


ae 


(11) 


A  typical  value  of  a/ is  150,  and  so  values  of  e  greater  than  0.2  correspond  to 


values  of  z  greater  than  about  30.  Equation  (9)  specifies  p  in  terms  of  E  . 

4 . 2  Changes  in  a  and  x 

The  changes  in  a  and  x  per  revolution,  denoted  by  Aa  and  Ax  respectively, 
are  given  by  equations  (12)  and  (13)  of  Part  I  as 


Aa 


-  2a26 


f * 

l  (1 


+  e  cos  E) 2 


e  cos  E) 


pdE 


(12) 


Ax  =  -  2a  6 


T.  } 

I  (cos  E  +  e)odE  , 

!  \  1  -  e  cos  E / 


(13) 


on  replacing  the  integrals  from  0  to  2tt  by  twice  the  integrals  from  0  to  it  .  Here 
6  is  a  constant  given  by  6  =  FSC^/m  ,  where  m  is  the  mass  of  the  satellite  and  F  is 
a  factor  which  allows  for  the  rotation  of  the  atmosphere.  From  Part  I, 


/,  rPow  V 

V  vd0  ' 


(14) 


pO 


where  v  is  the  velocity  of  the  satellite  relative  to  the  Earth's  centre  and  i  the 
inclination  of  the  orbit  to  the  equator. 


EVALUATION  OF  Aa  and  Ax 


When  e  >  0.2  ,  the  integration  of  (12)  and  (13)  can  best  be  performed  by  express¬ 
ing  E  in  terms  of  \  ,  where 
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E  =  1  -  A  /z  , 


so  that 


sin  E  =  A 


(§)'(’  -  trf 


JE  ■  (t)'('  - 11)  d>  • 


Substituting  (15)  into  (9)  gives 


)  bx2AA 

p  =  Ppjl*  — 


|  exp(-  A2)  , 


which  may  be  expressed  in  terms  of  p  by  using  equation  (5),  as 


[l  +  ipA4 1 1  -  JuA2(a2  -  2)|  j  exp(-  A2)  . 


We  now  make  the  simplifying  assumption  that  the  terms  of  order  p  in  (19)  can  be 
ignored,  and  take  as  our  standard  equation  for  density 

P  =  Pp(l  +  ihAA)  exp(-  A2)  .  (20) 

From  (19)  the  error  in  equation  (20)  is  hp  p2\6(\2  -  2)exp(-  A2)  .  For  heights  up 
l  2  \  P  '  2  ' 

to  2.4  H  above  perigee  l  A  <  2.4)  this  is  less  than  Q.125p  p  ,  so  that  the  error  in  (20) 
p  v  /  p 

is  less  than  0.005p  for  p  <  0.2  .  At  heights  more  than  2.4  H  above  perigee,  this 
P  P 

error  can  be  larger,  having  a  maximum  of  0.025p  at  4.7  H  above  perigee  for  p  =0.2  : 

P  P 

at  such  great  heights,  however,  the  assumption  of  constant  p  has  itself  become  invalid, 
and  it  is  pointless  to  pursue  the  effect  of  O^j2)  terms  into  these  rarefied  realms.  To 
summarize,  we  are  seeking  a  good  approximation  to  the  density  in  the  region  from  perigee 
up  to  about  2.4  H  above,  without  being  too  much  concerned  about  greater  heights,  where 

P  2 

the  error  could  increase  to  0.6p  p  ,  if  p  remained  constant. 

P 

If  we  substitute  for  cos  E  by  (15),  dE  by  (17)  and  p  by  (20),  we  find  that 

2 

equations  (12)  and  (13)  can  be  expanded  in  powers  of  A  /z  to  give 


-  2a  p  6 
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,  3 

+  e)?  [ 


2  -  l)*2 


i8e  -  3e  -  1 
4z ( 1  -  e2) 


0(4)U  ♦  ijjAA)exp(-  A2)dA 
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e  specified  in  Part  III.  Since  /2z  >  6  for 
(22)  may  be  replaced  by  “  ,  with  relative  error 


exp^-  X2)dX  =  , 


2  * 
a  p  6 
P 


_  3  +  e-  3  15mU  ♦  Q  +  ±\ 

UJTTT^lj1  8z(,  -e2)  ^  64.(,  -e2)  \Vf 


2  2 

Here  two  terms  inside  the  curly  brackets,  3K1/4z  in  (25)  and  3K2/4z  in  (26),  have 

been  omitted  because  it  was  found  in  Part  III  that  they  are  less  than  0.00005  if 

H/r  =  0.008  ,  for  0.2  e  <  1  . 

P 

6  PERIGEE  DISTANCE  AS  A  FUNCTION  OF  e 

In  Part  III  it  was  shown  that  when  equation  (25)  with  *j  =  0  is  divided  by 
equation  (26)  with  m  =  0  ,  the  0^1/z2)  terms  are  less  than  0.15/z2,  that  is,  less  than 
0.00016.  Hence 


4a  _  1  ~  e  3m ( 1 

Ax  2  z ( 1  +  e )  4z(1 


-  e)  ^  JO. 35m2  0.15\ 

+  e)  ^  z  ~z/ 


where  the  0^0.15/z2)  term  is  taken  as  incorporating  a  term  O^m/z2)  >  and  0.35m2/z  is 
the  maximum  value  of  this  term  (when  e  =  0.2  and  z  =  30  ) .  If  we  write 

H,  -  Hp(t.|„)  ,  (28) 

so  that,  as  (4)  shows,  Hj  is  the  value  of  H  at  a  height  f  H^  above  perigee, 
equation  (27)  reduces  to 
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Thus,  as  with  eccentricities  between  about  0.03  and  0.2,  the  constant-H  equation 

for  da/dx  may  be  used  if  H  is  evaluated  at  a  height  §  above  perigee.  Conversely, 

if  an  orbit  is  accurately  determined  from  observations  and  accurate  values  of  da/dx,  e 

and  x  are  available,  equation  (29)  can  be  used  to  evaluate  H  at  a  height  §  H  above 

perigee . 

This  is  a  neat  and  useful  one-revolution  result;  but,  in  order  to  integrate  (29) 

over  the  whole  lifetime  of  the  satellite,  we  have  to  take  account  of  the  fact  that 

is  evaluated  at  a  height,  y^  say;  which  is  slowly  decreasing  as  r^  decreases.  We 

need  a  constant  'mean'  value  for  and  fortunately  an  appropriate  value  can  easily  be 

guessed.  If  we  assume  that  e^  =  0.8  is  the  highest  practicable  value  of  e  ,  the 

maximum  decrease  in  perigee  height  is,  from  Fig  3  of  Part  III,  0.5H  c;  0.5H^g  (when 

e^  =  0.8  and  e  =  0.2  ).  As  the  orbit  spends  more  time  at  high  eccentricity  than  at  low 

(see  Fig  8  of  Part  III),  the  mean  perigee  height  is  likely  to  be  about  0.2H^^  below  the 

initial  perigee  when  the  orbit  is  contracting  from  e^  =  0.8  to  e  =  0.2  .  At  the  other 

extreme,  when  the  orbit  contracts  very  little,  the  mean  perigee  height  is  virtually  equal 

to  the  initial  perigee  height.  Thus  the  mean  decrease  in  perigee  height  during  the  life, 

for  all  varieties  of  orbit,  is  likely  to  be  between  0  and  0.2  H  y  and  may  be  taken 

as  0.1  H  „  with  maximum  error  0.1  H  n  .  Hence  the  most  promising  mean  height  for 
pu  pu 

evaluating  H  is  at  a  height  1.4  H  above  the  initial  perigee.  So  we  define 

p0 

H*  =  H  .(1  +  1.4jj)  ,  (30) 

pO 


and  try  it  as  a  mean  value  for  . 
Dividing  (28)  by  (30),  we  find 


H*  1  +  p(0.1  -  a)  -  p2(1.5o  +  0.14)  +  0 


(p3)j  * 


where 


rPo  ~  rp 
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The  variation  of  o  with  e/e^  as  given  by  Part  III  is  shown  in  Fig  2:  we  need  a  simple 
approximation  for  ci  in  terms  of  e  and  we  take 


■  •  («•*  -  -  0 


±  0.065 


which  is  shown  by  the  broken  lines  in  Fig  2.  The  error  quoted  is  the  maximum  for 
0.2  <  <  0.8  ,  and  occurs  when  =  0.8  and  e/e^  =  0.4  .  In  view  of  (31)  and  (33), 

equation  (29)  may  be  written 


V  v  V 

Oa.V.'.n'] 


*.  .■  , 


v.vn 


to 


=  1  +  _ 
dx  1  2( 


?[•  *  «)«-'  -  (■>•«  -  °-3eo0  -  ,)j]*.^  .  2^) 


(34) 


as  the  p  term  in  (31)  and  the  error  term  in  (33)  lead  to  0-terms  smaller  than 
0^0.35)j2/zj  .  In  (34),  we  need  to  express  e  in  terms  of  x  .  Since  the  perigee 
distance  decreases  by  less  than  jH*  , 


a  =  Qq  +  x  +  0($H*)  , 


(35) 


where  denotes  the  initial  perigee  distance, 


From  (35), 


Q0  =  ao(’  ■  eo)  =  a0  ■  xc 


(36) 


I  ■  • 


(37) 


Eliminating  e  from  (34)  with  the  aid  of  (37),  we  find 

1  d(a  -  x)  =  [l  *  *  (°-4  '  Q'^oX1  ~  e0)}_ 


H*  dx 


2(Qq  +  2x)x 


Qq  ^0Qq(0-4  -  °-3eo) 

2(q0  ♦  2x)x2 


(38) 


Integrating  (38)  and  collecting  terms,  we  have 


_E° _ P  . 

H* 


..  ‘oh,*2*), 

iln  r?73 — r  -  v  +  u 


+  2x0 ) 


(o . 5  +  O.Je  -  0.3e2)i{ln  -?1— - T 

V  0  on  ?(Q0  +  2x0j 


-  -  0-3eo)(J-ra  - 


+  0l0.35)j‘'ln 


0  0.15 


(39) 


The  term  in  curly  brackets  may  be  written  as  a  +  0(p)  ,  which  in  turn  may  be  expressed 
in  the  form  (32),  giving  finally 


pv;  p  _  i ,  __  f  v  \  u  j 
h*  "  u"irrr7TT- 


*  ttH  -  w(‘  -  *«)l5*o  -  'Af  -  <)  * »("-  •  tt)  ■ 


where  the  O^q2^  term  appears  because  ln^x^/x j  can  be  as  large  as  2.77  (when  e^  =  0.8 
and  e  =  0.2  )  and  2.77  x  0.35  =  0.97  .  The  coefficient  of  q  in  equation  (40)  has  a 
maximum  value,  when  e  =  0.2  ,  of  O.Ole^  (4  -  3Cq)(3Cq  -  l)(5eQ  -  1)  ,  which  is  always  less 
than  0.054  if  0.2  <  e^  <  0.8  ,  and  is  less  than  0.022  for  0.2  <  e^  <  0.6  .  The  q 
term  in  (40)  is  thus  of  order  0.05q  =  0.25q2  if  q  =  0.2  ,  and  can  be  absorbed  in  the 


o(„2) 


Equation  (40),  shorn  of  its  q  term,  shows  that  the  zero-jj  theory  can  still  be 
used,  with  error  0^q~^  »  if  H  is  evaluated  at  a  height  1.4  H  ^  above  the  initial 
perigee  height. 

Note  also  that,  from  equation  (23)  of  Part  III,  equation  (40)  may  be  rewritten  as 


IpO _ tp 


=  i  In 


;o(1  +  e)  , 

0/H  u2  0.15\ 

:(  1  +  eA\ 

°Ur  ’  p  ’  z  J  • 

[  0) 

\  p  / 

(=:  0.002) 

is  normally  less  than  q2 

VARIATION  OF  e  WITH  t 
From  (25)  and  (26),  we  have 


Ax  -  eAa  =  -  a  p  6 

P 


eiy 

\  ae  / 


(1  +  e):(1  -  e) 


!  ,  3  +  4e  -  3e2  .  3 

1 - 7 - 5T-  +  F  M  +  0, 

8z(l  -  e2  8 


(*  '  ?)]  ' 


The  1/z  term  in  (42)  has  the  value  (0.34  ±  0.04)y|z(1  -  e)|  for  0.2  <  e  <  0.8  , 
which  is  of  order  0.01  or  less  for  z(1  -  e)  >  30  .  The  term  is  dropped  after  integra¬ 
tion  in  Part  III,  and  here  it  is  more  convenient  to  drop  it  at  once.  On  dividing  by 


T  =  T, 


>(a/ao)1  * 


where  T  is  the  orbital  period,  equation  (42)  gives 

de  i  U  \!  .jiVUl  +  erO  ,  3  .  ^0.3  p  I  \1 

3T  -  t^(2”o)  pHp\T/ - 2 - ['  *  8  “  *  °(~  •  I  •  7)J  •  1 

To  proceed  further,  p  H^(a_/a)  must  be  expressed  as  a  function  of  e  .  On 

p  p  \  0  / 

taking  the  reference  distance  r^  in  equation  (3)  as  the  initial  value  r  ^  ,  we  may 

express  c  in  terms  of  o  in  the  form 

p  po 


Ppoj1  +  ba2H*2|exp(oH*/Hp0)  , 


(45) 


where  a  is  given  by  (32).  In  the  subsequent  analysis  the  O^p2^  terms  will  in  the 
end  be  neglected;  but  there  are  many  of  them,  and  we  need  to  evaluate  them  numerically 
to  ensure  that  they  do  not  accumulate  to  give  an  unacceptably  large  0-term.  The  fuller 
equation  for  b  in  terms  of  p  and  a  is  therefore  needed,  and  from  (5)  it  is 


b  =  ^2  j1  ‘  P(°  "  ^  +  °(p2)(  • 


(46) 


On  substituting  for  b  from  (46),  eliminating  H*/Hpg  by  means  of  (30)  and  expanding 
exp(1.4pa)  as  a  power  series,  we  can  rewrite  equation  (45)  as 


Pp  pp0 


jl  +  pa(1. 4  +  0.5a)  +  p2a2^2.4  +  0.2a  +  0.25a2j  +  O^p^jexp  a  .  (47) 


From  (6)  and  (30)  we  find 


K*  =  H*  jl  -  0.5pa  -  p2a( 0.7  +  0.125o)  + 
P  po 


°d)|  ■ 


(48) 


Multiplying  (47)  by  (48)  gives 


P  H1  =  p  n 
P  P  P° 


RJo| 


1  +  pa(0.9  +  0.5a)  +  0 


(°V)j 


exp  o  , 


(49) 


or,  on  using  (41)  to  evaluate  exp  a  , 


p  h*  =  P„n 

p  p  p° 


Hpo)I°(l'.V0)j  '  *  “{’  *  °-5<«  -  °-2>}  *  °(“2  •  JTJ  •  nr)  ■  <50> 


Here  0  <  a  <  0.5  ,  and  so  |0.5a(a  -  0.5) |  <  0.075  ,  giving  a  term  of 

0(0.075p)  =  0^0.375(j2j  if  p  =  0.2  :  this  term  can  be  absorbed  in  the  0^p2J  term  in 

(50).  Also,  from  (35) 


_£ 

a 


1  -  e 


1  -  e 


0 


<k)- 


(51) 


On  substituting  (50),  without  the  small  term,  and  (51)  into  (44),  and  using  (33)  to 
approximate  for  a  ,  we  obtain 


+  This  corrects  an  error  in  the  formula  for  b  in  Part  IV,  where  the  a  and  a  are 
interchanged . 


no 

vs 


is.  ,  .  _.V6  f'YoVVc  *  »)2C  --S)jji  ,  „(0.4  -  0.3e„)(- 

dt  TqO  -  e0)\  1  +  eQ  /  e  |  8  0/\ « 

*  42  •  ¥)j  • 

after  dropping  0-terms  that  are  smaller  than  others. 

Equation  (52)  now  has  to  be  integrated.  As  in  Part  III  we  write 

_  “OpO8  gv o^o\‘ 

'  To('  -  eo)  '  ’  *  »o  /  ’ 


so  that  (52)  becomes 


where  J  and  K  are  constants  given  by 


J  =  0.3eQ  -  0.025 


K  =  (o.4  -  0,3eQ)< 


Preparing  for  integration,  we  rewrite  (54)  as 

-£i|<,.Ju)j,.o(„2,¥)!  ■ 


(1  +  e)  (1  -  e)2  (1  +  er(1  -  e)' 


On  integration,  (57)  gives 
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Equation  (58)  with  p  =  0  should  be  the  same  as  equation  (4.126)  of  Ref  11.  Making  the 
comparison,  we  see  that  the  term  0(0. 3/z)  in  (58),  which  is  evaluated  numerically  in 
Part  III,  can  be  replaced  by  0^H/5r^j  =  0(0.002)  . 
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with  t  ,  but  in  order  to  bring  it  to  the 


t^  ,  the  lifetime  of  the  satellite.  In 


Equation  (58)  gives  the  variation  of  e 
same  format  as  in  Part  III,  we  need  to  define 
Part  III,  t  was  defined  as  the  value  of  t 

L 

t  in  terms  of  e  was  not  valid  for  e  <  0.2 
satisfactory  for  zero  jj  ,  but  causes  difficulties  when  jj  is  non-zero  because  the  K(j 
term  in  (57)  becomes  large  relative  to  the  'main'  term  as  e  -*■  0  .  As  the  definition  of 


for  e  =  0  ,  even  though  the  equation  for 
.  This  arbitrary  definition  was  quite 


t  is  arbitrary,  and  it  needs  to  be  approximately  the  satellite's  lifetime,  it  should  be 

L 


defined  in  such  a  way  as  to  be  consistent  with  the  theory  for  e  <  0.2  .  This  aim  can  be 
achieved  by  defining  t^  as  the  value  of  t  when  e  =  0  in  equation  (58)  with  the  last 
of  the  three  terms  in  Kjj  omitted.  We  then  have 


3  + 


(1  +  Jp)CtT 


(’  +  eo)i/(1  "  eo ) 
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F71] 

'1 

l1  +  e. 
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Kp 


1  +  3e,. 


v +  eoyc  -  eo) 


- 1 


(59) 


Subtracting  (58)  from  (59)  gives: 
(1  +  Jp)C 


tL  -  t  lj  +  O^j2  ,  =  ~  KpjjUe)  +  L^eQJ  -  L(e) 


(60) 


where 


f(e)  = 


3  +  e 


(1  +  e) /( 1  -  e) 


3 _ 1_  /2  +  /(I  -  e) 


/2  (/2  +  1)/(1  -  e) 


(61) 


as  in  Part  III,  and  L(e)  and  M(e)  are  new  functions  defined  by 

/2  +  /(I  -  e) 


L(e) 


In 


2/2  /2  -  /(I  -  e)  ’ 


(62) 


M(e)  = 


1  +  3e 


(1  +  e) /( 1  -  e) 


-  1 


(63) 


As  a  guide  to  orders  of  magnitude,  some  values  of  f(e),  L(e)  and  M(e)  are  listed  in 
Table  1. 


Table  1 


Values  of  f(e),  L(e)  and  M(e) 


e 

0.2 

0.4 

0.6 

0.8 

f  (e) 

0.0775 

0.3235 

0.8405 

2.112 

L(e) 

1.581 

1 .305 

1 .021 

0.6946 

M(e) 

0.4907 

1 .029 

1 .767 

3.224 
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r. 


Putting  e  =  and  t  =  0  in  (60)  gives 


(1  ♦  J^)CtL  =  f(eQ)  -  KyM(e0)  . 


(64) 


Dividing  (60)  by  (64),  we  have 


1  -4- 


f  (e) 

"  Kp 

|M(e)  -  L(e) ] 

[  KpL| 

w 

1 

fl 

E^y 

-  KpM(e0] 

1  £(V 

I  " 

KpM(eQ) 

(65) 


Equation  (65)  gives  t/t  in  terms  of  e  and  e  ,  and  Fig  3  shows  how  e/e 


varies  with  t/t 


0 

for  =  0.4,  0.6  and  0.8  with  p  =  0  and  p  =  0.2  .  The  diagram 


shows  that  the  effect  of  p 


L 

3.' 

on  t/t 


is  small  and  usually  negligible.  However,  this 


does  not  mean  that  p  has  a  negligible  effect  on  e  :  equation  (52)  shows  that  de/dt 
is  increased,  initially,  by  a  factor  (1  +  3p/8)  -  that  is  1.075  if  p  =  0.2  .  If  this 
factor  stayed  constant,  t  would  be  reduced  by  a  factor  1.075  while  the  variation  of  e 

Li 

with  t/t  would  be  independent  of  p  ;  in  reality,  however,  the  factor  slowly  increases 


as  time  goes  on,  thus  producing  the  small  variation  of  e/e^  with  p 
8  LIFETIME  IN  TERMS  OF  THE  INITIAL  DECAY  RATE 


seen  in  Fig  3. 


To  obtain  the  lifetime  t^  in  terms  of  T  we  go  back  to  equation  (25),  but  drop 
the  1/z  term  within  curly  brackets,  which  is  of  0^H/2r^  ,  and  the  p/z  term,  which 
is  smaller.  Then  we  have 
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(66) 


This  may  be  written  with  the  aid  of  (53)  as 

3CT 


(67) 


Multiplying  (67)  by  (64),  we  find 
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(68) 


where  0.2p  represents  the  maximum  value  of  the  neglected  cross-multiplied  p  terms. 
This  is  the  required  equation  for  the  lifetime  in  terms  of  the  initial  decay  rate. 


We  first  need  to  check  that  it  is  consistent  with  the  low-eccentricity  theory, 
the  theory  of  Part  I\ 
eccentricity  thee  y  gives 


From  the  theory  of  Part  IV,  as  revised  in  Ref  12,  we  find  that,  when  e^  =  0.2  ,  the  low- 


(low-e  theory)  . 
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u  o 


1.25 


1  -  0 . 1 8  7jj  +  0 


(69) 


When  equation  (68)  is  evaluated  at  e^  =  0.2  ,  we  find 
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Vo 

4Trt 


1.248 


1  -  0.091m  + 


°(^ 
'  D 


(70) 


2 

These  differ  by  0.096p  =  0.5jj  if  (j  =  0.2  and  are  therefore  consistent.  Thus 
equation  (68)  is  acceptable  even  for  e^  =  0.2  ,  when  the  whole  lifetime  is  outside  the 
validity  of  the  high-eccentricity  theory.  It  is  a  fortiori  acceptable  for  higher  values 
of  ,  because  the  added  time  from  e  =  e^  down  to  e  =  0.2  will  be  given  correctly 

by  the  high-e  theory. 


On  evaluation,  the  term  in  curly  brackets  in  (68)  proves  to  have  a  value  between 
0.035  and  0.091  for  0.2  £  e^  ^  0.8  .  Thus  we  can  say  that  the  constant-H  lifetime 
formula,  with  p  =  0  in  (68),  still  applies,  to  0^0.5p^j  . 

9  DISCUSSION  AND  CONCLUSIONS 


The  analysis  of  section  6  shows  that  the  constant-H  theory  for  the  variation  of 
perigee  height  with  eccentricity  can  still  be  used  when  H  varies  with  height,  if  H  is 


evaluated  at  a  height  1.5  above  the  perigee  height  y^  (for  single-revolution 
results)  or  1.4  H  „  above  initial  perigee  height  y  n  (for  long-term  results).  With 

pU  pU 


e/eQ  ,  and  Fig  4 


this  refinement.  Fig  3  of  Part  III,  giving  perigee  height  in  terms  of 
of  Part  III,  giving  T/T^  against  -  ej  ,  will  remain  valid.  These  results  show 


that,  with  high-eccentricity  orbits,  simple  exponential  density  models  can  still  be 
successfully  used  if  the  choice  of  value  of  scale  height  is  adjusted  appropriately.  In 
the  analysis  the  value  of  p(=  dH/dy)  is  assumed  to  be  less  than  0.2,  corresponding  to 


perigee  heights  y  >  160  km  for  low  solar  activity  or  y  >  190  km  for  high  solar 
P  P 


activity . 


The  analysis  of  section  7  specifies  the  effect  of  non-zero  p  on  the  variation  of 
eccentricity  with  time.  Equation  (52)  shows  that  initially  the  rate  of  decrease  of  e 
increases  by  a  factor  of  ( 1  +  0.375p).  This  factor  becomes  gradually  larger  as  time  goes 
on  and  at  its  maximum,  for  an  orbit  with  e^  =  0.77  and  e  =  0.2  ,  reaches  (1  +  0.857p) 
When  time  t  is  expressed  as  a  fraction  of  the  total  lifetime  t  ,  and  e  is  plotted 

L 

against  t/t^  >  c^e  effect  of  p  is  generally  negligible  (see  Fig  3).  Consequently  the 
variation  of  T/Tq  with  t/t^  is  also  virtually  the  same  as  in  Fig  10  of  Part  III. 

The  curves  of  e/e^  versus  t/t^  exhibit  quite  strong  curvature,  and  an  orbital 
parameter  having  a  nearly  linear  variation  with  t  would  be  welcome.  As  it  happens, 
z/Zq  fulfils  this  requirement  well,  and  Fig  4  shows  how  z/z^  varies  with  t/t^  • 

The  full  formula  for  lifetime  in  terms  of  the  initial  decay  rate  is  given  by 
equation  (68).  On  evaluating  the  p  terms  in  this  equation,  it  is  found  that  they  are 


of  order  0.5p  .  Thus,  although  it  is  more  accurate  to  use  equation  (68),  the  zero-p 

2 

equation  for  lifetime  is  essentially  valid  to  order  p 

In  regard  to  past  work  the  main  implication  of  this  Report  is  in  confirming  ideas, 
previously  only  surmised,  (a)  that  H  should  be  evaluated  at  a  height  |H  above  perigee 
in  orbit  determination,  and  (b)  that,  when  scale  height  is  determined  from  orbit  analysis 
over  a  short  time  interval,  this  is  the  height  at  which  H  applies. 

The  chief  value  of  the  theory  for  future  work  is  in  allowing  more  secure  and 
effective  use  of  simple  atmospheric  models  in  programs  for  orbit  determination  and 
analysis.  The  Report  shows  how  best  to  interpret  these  models  in  order  to  achieve  more 
accurate  numerical  results. 
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Fig  1  The  variation  with  height  y  of  density  scale  height  H  and  its  gradient 
M  =  dH/dy  as  given  by  CIRA  1972.  for  high  and  low  solar  activity 


